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Abstract 

A non-classical Weyl theory is developed for skew-self-adjoint Dirac 
systems with rectangular matrix potentials. The notion of the Weyl 
function is introduced and direct and inverse problems are solved. 
A Borg-Marchenko type uniqueness result and the evolution of the 
Weyl function for the corresponding focusing nonlinear Schrodinger 
equation are also derived. 
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1 Introduction 

Consider a skew-self-adjoint Dirac-type (also called a Dirac, ZS or AKNS) 
system 



A. 

dx 



y{x, z) = {izj + jV{x))y{x, z) (x > 0), 
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where 1^^. is the x ruk identity matrix and v{x) is an mi x m2 matrix 
function, which is called the potential of the system. Skew- self-adjoint Dirac- 
type systems are very well known in mathematics and applications, though 
the Weyl theory of such systems is non-classical. The name ZS-AKNS comes 
from the fact that system (1.1) is an auxiliary linear system for many im- 
portant nonlinear integrable wave equations and, as such, was studied, for 
instance, in [1,3,8,11,14,18,29,33,42,48]. For the case that mi ^ m2 systems 
of the form (1.1), (1.2) are, in particular, auxiliary linear systems for the cou- 
pled, multicomponent, and mi x m2 matrix nonlinear Schrodinger equations. 
We note that the above-mentioned case mi 7^ m2 was much less studied. In 
this respect our present paper is a continuation of the research started in the 
papers [12, 13], where self-adjoint Dirac systems with rectangular potentials 
were discussed. We solve direct and inverse problems, that is, we construct 
Weyl functions and recover the mi x m2 potential v from the Weyl function, 
respectively. 

Spectral and Weyl theories of auxihary linear systems form a basis of the 
Inverse Spectral Transform (ISpT) approach [5,6,19,29,30,32,35,38,40,41] 
to several important initial-boundary value problems for integrable nonlinear 
equations. In particular, the evolution of the Weyl functions in the form of 
linear fractional transformations (see (4.15)) was first proved in [38,40]. The 
case of the focusing matrix nonlinear Schrodinger (fNLS ) equation 

"^vt + i(v^x + '2vv*v) - 0, 

where v are square matrices, was studied in [29], and here we derive the 
evolution of the Weyl functions for the fNLS with rectangular matrices. 

As usual, N stands for the set of natural numbers, M stands for the real 
axis, M_|_ stands for the positive real axis, C stands for the complex plain, 
and C_|_ for the open upper semi-plane. The notation Cjv^ stands for the open 
semi-plane {z : ^{z) > M > 0}, where ^{z) (or ^z) is the imaginary part of 
z. If a e C, then $R(a) is its real part and a its complex conjugate. If £ is a 
set, then C is the closure of the set. The notation Im is used for image. An 
m2 X mi matrix a is said to be non-expansive, if a*a < Imi (or, equivalently, 
if aa* < /ma)- We assume that the spaces C™' are equipped with the P- 
norm and use ||q!|| to denote the corresponding operator norm of the matrix 
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a. If the matrix function / is absolutely continuous and sup ||/'|| < oo, we 
say that / is a function with bounded derivative or (equivalently) that / is 
boundedly differentiable. The set of functions boundedly differentiable on 
[0, /] is denoted by -B^[0, l\. We use / to denote the identity operator and 
i?(Hi, H2) to denote the class of bounded operators acting from Hi to H2. 
We write 5 (Hi) if Hi = H2. 

We set mi + —: m. The fundamental solution of system (1.1) is 
denoted by u{x, z), and this solution is normalized by the condition 

U{0,z)^lm. (1.3) 

2 Direct problem 

We consider the Dirac system (1.1) on the semi-axis x G [0, 00) and assume 
that V is measurable and locally bounded, that is, bounded on all finite 
intervals. Using a method which is similar to the method proposed in [29,33], 
where the case mi = 1712 was treated, we shall apply Mobius transformations 
and matrix balls to solve the direct problem for Dirac system. 

We next introduce a class P>(j, Cm) of nonsingular m x mi matrix func- 
tions V(z) with property-j, which are an immediate analog of the classical 
pairs of parameter matrix functions. More specifically, the matrix functions 
V{z) are meromorphic in C_|_ and satisfy (excluding, possibly, a discrete set 
of points) the following relations 

V{zyV{z) > 0, V{z)*jV{z) > ioT z e Cm ^ {z : "^z > M > Q}. 

(2.1) 

Definition 2.1 The setAf{x,z) of Mobius transformations is the set of val- 
ues {at the fixed points x, z) of matrix functions 

ip{x,z,V)^[Q Im,]u{x,z)-^V{z)(^[lm, o] ii(x, z)-^7'(^)) (2.2) 
where V{z) are nonsingular matrix functions with property-j, that is, 

V{z)eP>{j,CM). (2.3) 
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Proposition 2.2 Let the Dirac system (1.1) on [0, oo) be given and assume 
that 1 1 I'll is bounded by M: 

\\v{x)\\<M fora;eM+. (2.4) 

Then the sets M{x,z) are well-defined in Cm- There is a unique matrix 
function ip{z) such that 

ifiiz)^ fl Af{x,z). (2.5) 

x<oo 

This function is analytic and non-expansive in Cm- 

Proof. The proof is similar to the proof of [12, Proposition 2.2], only z is 
considered now in the semi-plane Cm instead of the semi-plane C+ in [12]. 
We shall not, however, skip the details of this important proof, although 
we will later omit the details of proofs similar to those in [12, 13]. It is 
immediately apparent from (1.1) and (2.4) that 

-^(u{x, z)*ju{x, z)) — —2u{x, zy{'is{z)Im — V{x))u{x, z) < 0, z e Cm- 
ax 

(2.6) 

According to (1.3) and (2.6) we have 

2l(x, z) = {^ij{x, z)}lj^-^ := u{x, z)*ju{x, z) < j, z E Cm, (2.7) 

where 21 is partitioned into four blocks so that is an rrii x nii matrix 
function {i — 1,2). Inequality (2.7) yields 

{u{x,zy)-'Mx,z)-'>j. (2.8) 

In view of (2.1), we note that V{z) maps C"*i into a maximal j-nonnegative 
subspace of C"*. Inequality (2.8) implies that u{x,z)~^ maps all maximal j- 
nonnegative subspaces into maximal j-nonnegative subspaces, that is, the 
image Im (^u{x, z)~^V{z)^ is a maximal j-nonnegative subspace. Clearly 
Im ( [/^^ O] ) is a maximal j-nonnegative subspace too. Thus, we get 

det([V 0]u{x,zy'V{z)^ ^0, (2.9) 
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and so M is well-defined via (2.2). Indeed, if (2.9) does not hold, there is a 
vector / e C"*i such that 

[V 0] ju{x, z)-'V{z)f = [I^, 0] u{x, z)-'V{z)f = 0, fy^O. (2.10) 

Since the left-hand side in (2.10) equals zero and (as discussed above) we 
deal with the maximal j-nonnegative subspaces, we have u{x, z)~^V{z)f G 
Im ( [imi O] ). But then it follows from the second equality in (2.10) that 
/ = 0, which contradicts the inequality in (2.10). 
Next, we rewrite (2.2) in the equivalent form 

= u{x,z)-'r{z)(^[l^, 0]u{x,z)-'V{z)y\ (2.11) 

In view of (2.1), (2.11), and of the definition of 21 in (2.7), the formula 

ip{z)eM{x,z) (2.12) 

is equivalent to 

[Im, ^(^)*]21(X,Z) 

In the standard way, using formula (2.6) and the equivalence of (2.12) and 
(2.13), we get 

J\f{xi,z)GN'{x2,z) for Xi>X2. (2.14) 

Moreover, (2.13) at a; = means that 

J\f{0,z) = {ifi(z) : ifi{z)*ifi{z) < Im,}. (2.15) 

By Montel's theorem, formulas (2.14) and (2.15) imply the existence of an 
analytic and non-expansive matrix function <^{z) such that 

ip{z)e f]j\f{x,z). (2.16) 

x<oo 

Indeed, because of (2.14) and (2.15) wc sec that the set of functions ip{x, z, V) 
of the form (2.2) is uniformly bounded in Cm- So, Montel's theorem is 



-'mi 

ip{x,z,V)^ 



3{z) 



> 0. 



(2.13) 



5 



applicable and there is an analytic matrix function, which we denote by 
(fiooiz) and which is a uniform limit of some sequence 



ipoo{z) = lim ip{xi, z, Vi) {i e N, Xi t, lim Xi = oo) 



(2.17) 



on all the bounded and closed subsets of Cm- Since Xi t and equalities (2.11) 
and (2.14) hold, it follows that the matrix functions 



Vijiz) ■.= u{xi,z) 



-'mi 

ip{xj,z,Vj)^ 



ij > ^) 



satisfy relations (2.1). Therefore, using (2.17) we see that (2.1) holds for 

If) 



Vi,oo{z) := u{xi, z) 



mi 



which implies that we can substitute V — Vi^oo and x — Xi into (2.11) to get 

'^^{z)eN{xi,z). (2.18) 

Since (2.18) holds for all i G N, we see that (2.16) is true for ip{z) = ipoo{z). 

Now, let us show that jV is a matrix ball. It follows from (2.6) and (2.7) 
that 

-^21 < -2{Q{z) - M)u*u < 2{Q{z) - M)2l < 2{Q{z) - M)j. (2.19) 

Taking into account the inequalities above and the equality 2l22(0, z) = —Im2-i 
we obtain 



- 2l22(x, z)>{l + 2{^{z) - M)x)lm, (z e Cm). 



(2.20) 



Note also that (2.7) implies Qi{x,z) ^ > j for z e Cm (see [27]). Thus, we 
get 



(2.21) 



1 /9 

Since — 2I22 > 0, the square root T = ( — 2I22) is well-defined and we 
rewrite (2.13) in the form 



- 2112212^2121 - {(p*T - 2I12T-1) {T(p - T-i2t2i) > 0, 



where 2ti2 = ■ Equivalently, we have 

(p = piUpr - 21^2^2121, U*U < Im2, (2-22) 

Pi := = ( - 2122)"'^', Pr (2lii - 2li22l2-2'2l2i)'^'. (2.23) 

Here uj is an m2 x nii non-expansive matrix function, pi and pr are the left 
and right semi-radii of the Weyl disc. Since (2.12) is equivalent to (2.22), 
the sets A/'(x, z) (where the values of x and z are fixed) are matrix balls. 
According to (2.20), (2.21), and (2.23) the next formula holds: 

Pl{x,z)^0 {X ^ oo), Pr{x,z)<Imi- (2.24) 

Finally, relations (2.16), (2.22), and (2.24) imply (2.5). ■ 

Remcirk 2.3 It follows from the proof of Proposition 2.2 that the inequality 

\\v{x)\\<M for 0<x</<oo (2.25) 

implies that the embedded matrix balls jV(a;, z) are well-defined for z e Cm 
and X on the interval [0, /]. 

Remcirk 2.4 If(p{z) G A/'(/, z) for all z e Cm and, in addition, (p is analytic 
in Cm, then it admits representation (2.2), where x — I andV satisfies (2.3). 
This statement is immediately apparent, since we can choose 



V{z) ^u{l,z) 



The property- j of this V follows from the equivalence of (2.12) and (2.13). 

In view of Proposition 2.2 we can define the Weyl function of the skew-self- 
adjoint Dirac system similar to the cases of the canonical system [41] and 
self-adjoint Dirac system [12]. 

Definition 2.5 The Weyl-Titchmarsh {or simply Weyl) function of Dirac 
system (1.1), where 11^(2;) || is bounded on M+, is the function tp given by 
(2.5). 
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Remcirk 2.6 From Proposition 2.2 we see that under condition (2.4) the 
Weyl- Titchmarsh function always exists and that it is unique as well. 

Now we introduce a more general definition of tlie Weyl function. For 
this reason we need the following corollary. 

Corollary 2.7 Let the conditions of Proposition 2.2 hold. Then the Weyl 
function is the unique function, which satisfies the inequality 



mi 



(p[z) \ U[X,Z) U[X,Z) 



-'mi 



dx < oo. 



(2.26) 



Proof. Using the first inequality in (2.19) we derive 



/ [imi ^{zy] 

Jo 



u{x, zyu{x, z) 



dx < 



2{{Q{z) - M) 



[imi v{zy] 



(2.27) 



X (21(0, z) - 2l(r, ^)) 



-'mi 



Recall that (2.12) imphes (2.13), and so for (/? satisfying (2.5) we get 



[I^, ^(z)*]2l(x,. 



mi 

ip{z) 



> for all < X < oo. 



(2.28) 



Inequality (2.26) follows immediately from (2.27) and (2.28). 
Moreover, since u*u > —21, inequahty (2.20) yields 



[O Im2\u{x,zyu(x,z) 





Jm2_ 



dx > rl. 



m2 ■ 



(2.29) 



In view of (2.29), the function (p satisfying (2.26) is unique. 



Definition 2.8 The Weyl function of Dirac system (1.1) on [0, oo), where 
\\v\\ is locally summahle, is the function </? satisfying (2.26). 

From Corollary 2.7, we see that Definition 2.8 is equivalent to Definition 
2.5 for the case that the potentials are bounded, but it can also be used in 
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more general situations [8]. A definition of the form (2.26) is a more classical 
one and deals with solutions of (1.1) which belong to I>^(0, oo). Compare 
(2.26) with definitions of the Weyl-Titchmarsh or M- functions for discrete 
and continuous systems in [22,23,41,45,46] (see also references therein). For 
the case of system (1.1), where mi = 1712, Definition 2.8 was introduced and 
the corresponding Weyl function was studied in [8,18,29,30,33]. 

3 The inverse problem on a finite interval 

In this section we consider a procedure for solving the inverse problem, that 
is, for recovering the potential v (and, hence, the Dirac system) from the 
Weyl function. We note that it suffices to recover v almost everywhere on all 
the finite intervals [0, /]. Therefore, we introduce the corresponding definition 
of Weyl functions of the skew-self-adjoint Dirac system on an arbitrary fixed 
interval [0, I] and solve the inverse problem on this interval. 

Definition 3.1 Weyl-Titchmarsh {or simply Weyt) functions of Dirac sys- 
tem (1.1) on [0, /] which satisfies (2.25), are the functions of the form (2.2), 
where x = I and V satisfies (2.3), i.e., the analytic functions ip{z) e Af{l, z) 

Recall that in view of Remark 2.3 the set A/'(/, z) is well-defined. 
We assume that v satisfies (2.25) and let 

I3{x)=[l^, 0]u{x,0), ^ix)=[0 lrn,]u{x,0). (3.1) 

It follows from sup^.^.; \\v\\ < 00 and from (1.1) that 

sup ||7'(a;)|| < 00, 7' := 4^7- (3-2) 

x<l ClX 

Moreover, from (1.1)- (1.3) we get 

u{x,z)*u{x,z) = u{x, z)u{x,z)* = Im- (3.3) 
Therefore, (3.1) implies 

^r = V, 11*^ Ira,. /37* = 0. (3.4) 
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Next, as in the self-adjoint case [13], we need the similarity result below for 
the Volterra operator 

r F{x)G{t) -dt, Ke B{Ll^{Q, /)), (3.5) 

^0 

where F{x) is an mi x m matrix function, G{t) is an m x mi matrix function, 
and 

iF{x)Gix) = Im,. (3.6) 

Proposition 3.2 [36] Let F and G be boundedly differentiable and let (3.6) 
hold. Then we have 

K^EAE-\ A-.^-iT-dt, A,E,E-'eB{Ll^{0,l)), (3.7) 

^0 

where K is given by (3.5) and E is a triangular operator of the form 

{Ef){x) = p{x)f{x) + E{x, t)f{t)dt, = iF'Gp, det p(0) ^ 0. 

(3.8) 

Moreover, the operators E^^ map functions with bounded derivatives into 
functions with bounded derivatives. 

We set 

F = 7, G^-iT*, 77* = /m., l^B\QJl (3.9) 

where stands for the class of boundedly differentiable matrix functions. 
Clearly F and G in (3.9) satisfy the conditions of Proposition 3.2. 

We partition ^{x) into two blocks 7 = [71 72] , where 71, 72 are, respec- 
tively, m2 X mi and m2 x m2 matrix functions. We show that, without loss 
of generality, one can choose E satisfying condition 

^"S = /m„ (3.10) 

where E~^ is applied to 72 columnwise. 
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Proposition 3.3 Let K he given by (3.5), where F and G satisfy (3.9), 
and let E be a similarity operator from Proposition 3.2. Suppose Eq G 



Then the operator E :— EEq is another similarity operator from Proposition 
3.2, which satisfies the additional condition (3.10). 

Proof. The proof of the proposition is similar to the case mi = m2 (see, 
e.g., [36, pp. 103, 104]) and coincides with the proof of Proposition 2.2 in [13]. 
■ 

Remark 3.4 The kernels of the operators E^^, which are constructed in 
[36], as well as the kernels of the operators Eq^ from Proposition 3.3 are 
bounded. Therefore, without loss of generality, we always assume that the 
kernels of E^^ are bounded. 

The next lemma easily follows from Proposition 3.3 and will be used to 
construct a representation of the fundamental solution u. 

Lemma 3.5 Let 7 he an m2 x m matrix function, which satisfies the last 
two relations in (3.9), and set 



S:^E-\E*)\ n:=[$i$2], G S(C-^ L^^(0, Z)); (3.12) 
{^ig){x) = $i(x)^, $i(x) := (£;-Si)(x); ^29 = Im,9 = 9] (3.13) 



where E is constructed {for the given 7) in Proposition 3.3 and g is used to 
denote finite- dimensional vectors. Then A, S, and U form an S-node, that 
is {see [37, 39, 4IJ) , the operator identity 








(3.11) 



AS - SA* 



-inn* 



(3.14) 



holds. Furthermore, we have 



00 




(3.15) 



i=0 
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Proof. Because of (3.5), (3.7) and (3.9), we get 

EAE-^ - {E-yA*E* ^K-K*^ -i-f{x) f ^{t)* ■ dt. (3.16) 

Jo 

The second equality in (3.12) and formulas (3.10) and (3.13) lead us to the 
equalities 

= (E-S) {x)g, n*/ = /' (i^-S) {tyfm = f ^{ty {{E-'y/) {t)dt. 

Jo Jo 

(3.17) 

Using the first equality in (3.12) and relations (3.17) we can rewrite (3.16) 
as the operator identity (3.14). To prove (3.15) we will show that 

N N N 

J]lm (^{A*yS-^U^ D ^Im {S-^Am) = J] Im {Am). (3.18) 

i=0 i=0 i=0 

For this reason, we rewrite (3.14) as S'^A = A*S-^ - iS-mU*S~\ Hence, 
for Ni, N2 > we obtain 

Im (^{A*f"^'s-'A^m^ (3.19) 

N1+N2 

+ Im (^{A*yS-m^ Dim (^{A*f'S-'A'''+my 

i=0 

Using (3.19), we derive (3.18) by induction. In view of (3.18), it suffices to 
show that 



J]lm(A^n)=L^^(O,0, (3.20) 

i=0 

which, in turn, follows from 11 = [$i $2] , where $25' = 9 (see (3.13)). ■ 

Remcirk 3.6 Given an S-node (3.14), we introduce a transfer matrix func- 
tion in Lev Sakhnovich form {see [37, 39,41]) '■ 

WA{r, z) := Im - izWS;^{I - zArY^Prli, < r < I, (3.21) 

12 



where I is the identity operator; A^., Sr G B(^L'^^{0, r)), 

A := PrAP;, Sr := PrSP;-, (3.22) 

A is given by (3.7); the operators S and U are given by (3.12) and (3.13), 
and the operator Pr is an orthoprojector from L'^^{0, I) on L^^{0, r) such 
that 

{Prf)ix) = fix) iO<x<r), feLl^iOJ). (3.23) 

Since P^A — PrAP*Pr, it follows from (3.14) that the operators A^, Sr, and 
Pril also form an S-node, i.e., the operator identities 

ArSr - SrA; = -iPrllll* P^ (3.24) 

hold. 

Now, in a way similar to [28, 29] the fundamental solution w of the system 

-^w(x,z) = —iz'y(xY'y(x)w(x,z), w(Q,z) = Im (3.25) 
dx 

is constructed. 

Theorem 3.7 Let 7 be an m2 x m matrix function, which satisfies the last 
two relations in (3.9). Then, the fundamental solution w given by (3.25) 
admits the representation 

w{r,z) ^ WA{r,z), (3.26) 
where WA{r, z) is defined in Remark 3.6 {see (3.21)). 

Proof. The statement of the theorem follows from the Continual Factoriza- 
tion Theorem (see [41, p. 40]). More precisely, our statement follows from 
a corollary of the Continual Factorization Theorem, namely, from Theorem 
1.2 [41, p. 42]. Using Lemma 3.5 we easily check that the conditions of 
Theorem 1.2 [41, p. 42] are fulfilled. (We note that, though (3.15) is vahd, 
this condition is not essential for the factorization Theorem 1.2 [41, p. 42] 
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because S is invertible.) Theorem 1.2 [41] implies: if 11*5'^ ^P^II is boundedly 
differentiable, then we have 

-^WA{r,z) = -izH{r)wA{r,z), lim WA{r,z) = Im, (3.27) 
ar r-^+o 

H{r):=-^{U*S;'PrU), (3.28) 

where wa is given by (3.21). Taking into account that E^^ are lower trian- 
gular operators, we derive 

PrEP;Pr = PrE, (E"^) *P; = P;P, (E"^) ^P,*. (3.29) 

Hence, formulas (3.12) and (3.22) lead us to 

S;^ = E^Er, Er := PrEP;. (3.30) 

Therefore, using (3.17), we rewrite (3.28) as 

i7(r) = 7(r)*7(r). (3.31) 

Compare formulas (3.25) and (3.27), (3.31) to see that (3.26) holds. ■ 

Now, consider again the case of a Dirac system. Because of (3.1) and (3.4) 
we obtain 

u{x,0)'y{x)*'f{x)u{x,0)~^ — 

In view of (3.25) and (3.32), direct differentiation shows that the following 
corollary of Theorem 3.7 holds. 

Corollary 3.8 Letu{x,z) be the fundamental solution of a skew-self- adjoint 
Dirac system, where the potential v is bounded, and let 7 be given by (3.1). 
Then u{x, z) admits the representation 

u{x, z) = e^^^'u(x, 0)w(x, 2z), (3.33) 

where w has the form (3.26) and the S-node generating the transfer matrix 
function wa is recovered from 7 in Lemma 3. 5. 





Im, 



(3.32) 
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Remcirk 3.9 For the case that j is given by (3.1), it follows from (1.1) and 
(3.4) that 

y (a;)7(a;)* = - [v{xy O] u{x, 0)7(0;)* = -v{xy /3(x)-f{xy = 0. (3.34) 

Thus, from (3.9) and (3.34) we see that F'G = 0, and so p in (3.8) is a 
constant matrix. Therefore, since 72(0) = Im,^ and $2(2;) = Im2> equality 
(3.17) implies that p{x) = 7^2, and formula (3.8) can be rewritten in the 
form 

{Ef){x)^f{x)+ r E{x,t)f{t)dt, EeB{Ll^{0,l)). (3.35) 
^0 

Important publications by F. Gesztesy and B. Simon [16, 17, 43] gave 
rise to a whole series of interesting papers and results on the high energy 
asymptotics of Weyl functions and local Borg-Marchenko-type uniqueness 
theorems (see, e.g., [7,9,10,13,15,21,31,33] and references therein). Here we 
use the 5'-node scheme from [29] and generalize the high energy asymptotics 
result from [33] for the case that Dirac system (1.1) has a rectangular mi x m2 
potential v, where mi is not necessarily equal to m2- We first recall that 5" > 
and $i(x) = (£'~"^7i)(x) is boundedly differentiable in Lemma 3.5. Next, we 
note that since 71 (0) = and E has the form (3.35), the equality 

$1(0) = (3.36) 

holds. Therefore, using Theorem 3.1 from [13] and substituting (3.36) there, 
we get the statement below. 

Theorem 3.10 Let U — [$i $2] be constructed in Lemma 3.5. Then there 
is a unique solution S e B(^L^^{0,1)^ of the operator identity (3.14). This 
S is strictly positive {i.e., S > 0) and is defined by the equalities 

pi /'min{x,t) 

S^I+ / s{x,t) -dt, s{x,t):^ / $;(x-C)$'i(i-C)*c?C- (3.37) 

Jo Jo 

Now, we can apply the ^-node scheme to derive the high energy asymptotics. 

Theorem 3.11 Assume that cp is a Weyl function {i.e., (p admits repre- 
sentation {2.2), where x — I and V satisfies (2.3)), and the potential v of 
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the corresponding Dirac system (1.1) is bounded on [0, I] {i.e., (2.25) holds). 
Then {uniformly with respect to ^{z)) we have 

^{z) = 2iz e^''''^i{x)dx + O (2ze^'^y^/^^ , ^{z) oo. (3.38) 



Proof. Step 1. Similar to the self-adjoint case, we first show that 

I. 



{I - 2zA)-^U 



mi 

_^{z) 



< C/VQz for some C > 0. 



For this purpose, we consider the matrix function 

U{z) = [I^, ip{zy] {wa{1, 2zywA{l, 2z) - Q 



^ mi 

^{z) 



(3.39) 



(3.40) 



The proof that this function is bounded is different from the self-adjoint case 
(compare with [13]). More specifically, because of (3.3), (3.26) and (3.33) we 
have 



i^(^) = e"(^-^) [/^, ^{zY]u{l,zru{l,z) 



- Jmi - ip{z)*ip{z) 



We substitute (2.11) into the formula above and rewrite it as 

U{z) =e^'(^-^)(( [I^, 0] u{l,z)-^V{z)Vr{z)-'y'y {V{z)Vr{z)-y 
X {V{z)V,{z)-'){[l^, 0] u{l,z)-'V{z)V,{z)-')-' 



'mi 



ip{z)*ip{z), 



(3.41) 



where Vi is the upper mi x mi block of V. Clearly, inequalities (2.1) imply 
the invertibility of Vi{z) and boundedness of V{z)Vi{z)~^ . 

Step 2. To derive (3.39) we also need to examine the asymptotics of 
u{l,z) {z — )■ oo). This is achieved by considering the procedure for con- 
structing the transformation operator for Dirac system (sec also some more 
complicated cases in [20,24,36]). We will show that u admits an integral 
representation 



J — 



ii(x, z) = e'^^^ + / e'^'N{x,t)dt, sup ||A^(x, i)|| < oo {0<\t\<x<l). 

(3.42) 
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Indeed, it is easily checked that 

oo 

u{x,z) = y^^Ui{x,z), 

i=0 

uoix, z) := e"'^^ Viix, z) := 



(3.43) 



^i{x-t)zj ■ 



jV{t)ui^i{t,z)dt for i > 0, 



since the right-hand side of (3.43) satisfies both equation (1.1) and the nor- 
mahzation condition at a; = 0. Moreover, using induction we derive in a 

standard way that 

Ui(x,z)^ f ^^'Ni(x,t)dt, sup \\Ni(x, t)\\ <(2Myx'-^/(i-l)\, (3.44) 

J —X 

where i > and M is the same as in (2.25). We see that formulas (3.43) and 
(3.44) yield (3.42). 

Step 3. Because of (3.3) and (3.42) we get 



u{l, z) ^ = u{l, z)* 



—ilz 



e^"^! 



1Tl2. 



+ o(l) ) {zeC+, z^ oo). 

(3.45) 



In view of (2.2), (3.41) and (3.45), for any e > there are numbers M > 
M > and C > 0, such that for all z satisfying Q^;^ > M we have 

Mz)\\<e, \\U{z)\\<C. (3.46) 

Next, using (3.14) and (3.21) we easily obtain (see, e.g., [15,39]) 

wa{1, zywA{l, z)^Im + i(^ - z)U*{I - zA*)-^S,^\l - zA)-^U. (3.47) 

We substitute (3.47) into (3.40) to rewrite the second inequality in (3.46) in 
the form 



2i(z - z) [Im, <p{z)*] n*(7 - 2zA*)-^Si\l - 2zA)-^n 



-'mi 



^ Clmi- 



(3.48) 



Since 5" is invertible, positive and bounded, inequality (3.48) yields (3.39). 
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Step 4. We easily check directly (see also these formulas in the works 
on the case mi = that 

(/ - zA)-^f = f-iz / e'^^-''^'f{t)dt, (3.49) 

Jo 

mi -2zA)-^f = ! e^'^''-^^'f{x)dx, (3.50) 
Jo 

mi - 2zA)-'^2 = ^{e-'''' - (3.51) 

Now, apply — i$2 to the operator on the left-hand side of (3.39) and use 
(3.50) and (3.51) to get (uniformly with respect to "^{z)) the equality 



(3.52) 



2^' Jo 

Because of (3.52) (and the first inequality in (3.46)), we see that (3.38) holds. 
■ 

The integral representation below follows from the high-energy asymptotics 
of (p and is essential in interpolation and inverse problems. 

Corollary 3.12 Let ip he a Weyl function of a Dirac system on [0, /] which 
satisfies the condition (2.25) of Theorem 3.11. Then we have 

= ie-l.i.m._£e-^||±i|c , > M. (3.53) 

where l.i.m. stands for the entrywise limit in the norm of L^{0,21). 

Proof. According to Definition 3.1 and relations (2.14) and (2.15), the 
matrix function (f is non-expansive in Cm- Since (p is analytic and non- 
expansive in Cm, it admits (see, e.g., [26, Theorem V]) a representation 

POO 

ip{z) = 2iz / e'^''''^x)dx, z = ^ + iri, 77 > M > 0, (3.54) 
Jo 

where e~^''^$(a;) G Lm^xmii^^ °°)- Because of (3.38) and (3.54) we derive 

rl /"OO 

/ e2'(=^-'>($i(x) -$(x))dx= / e2^(^-')^$(x)dx + 0(l/v^(i)) 
Jo Ji 
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for ^{z) +00. Taking into account that e'^^^+^^^x) e L'^^y.rnii^, 00) 
for £ > 0, we rewrite the formula above as 

e2^(^-')"($i(a;) - ^x))dx = 0(1/7^), (3.55) 

Jo 

and the equahty (3.55) is uniform with respect to ^{z). Clearly, the left-hand 
side of (3.55) is bounded in the domains, where '<s{z) is bounded from above. 
Hence, in view of (3.55), its left-hand side is bounded in all of C and tends 
to zero on some rays. The identities 

[ e2'(^-')^($i(x) - ^x))dx = 0, i.e., ^i{x) = $(x) (3.56) 
Jo 

are clear. 

Using the Plancherel Theorem, we apply the inverse Fourier transform 
to formula (3.54) and get a representation of the form (3.53) for <l>. Since, 
according to (3.56), we have $1 = $, the formula (3.53) is valid. ■ 

Directly from (1.1) and (3.1), we easily get a useful formula 

P'ix) = v{x)-f{x), (3.57) 

which, because of (3.4), imphes 

/3'/3* = 0, /3V = V. (3.58) 

Thus, we will solve the inverse problem and recover v from (p once we recover 
/3 and 7. The recovery of $1 from cp is studied in Corollary 3.12. The next 
step to solve the inverse problem is to recover ^ from $1 in the following 
proposition. 

Proposition 3.13 Let Dirac system (1.1) on [0, I] satisfy the conditions of 
Theorem 3.11 and let 11 and S he the operators constructed in Lemma 3.5. 
Then the matrix function (3, which is defined in (3.1), satisfies the equality 

m-[lm, 0] - £ [S-'^[yty [^,{t) I,r^,]dt. (3.59) 
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Proof. First, we use (3.10) and (3.13) to express 7 in the form 

7(x) = (E[$i Im,]){x). (3.60) 

It follows that 

i{EAE-^E^[){x) = 71 (x) - 72(a;)$i(+0). (3.61) 

According to (3.1), (3.35), (3.60) and Remark 3.4, we have 

7i(0) = 0, $i(+0) = 0. (3.62) 

Therefore, we rewrite (3.61) as 

i{EAE-^E^[)(x) = 7i(x). (3.63) 

Next, we substitute K = EAE"^ from (3.7) into (3.63) and (using (3.5) and 
(3.9)) we get 

71 (x) = 7(^) £ lity {e^[) {t)dt. (3.64) 

Formulas (3.60) and (3.64) imply 

^iix)=^{x)J\E[^, /„,])(t)*(E$;)(t)cit. (3.65) 

Because of (3.29), (3.30) and (3.65) we see that 

^{x)e{xy = for e{x):^[lm, 0]-^ {E^[){ty (3.66) 

x(E[$i Im,]){t)dt=[lm, 0]-^ (5,i$;)(t)*[$,(t) i^^]dt, 

where := P^SP*. We will show that 9 ^ p. 

In view of (3.60) and the second relation in (3.66), we have 

e'{x) = - (a;)*7(x). (3.67) 

Therefore, (3.4) leads us to 

I3{x)e'{xy = 0. (3.68) 
20 



We furthermore compare (3.4) with the first relation in (3.66) to see that 

9{x) = h{x)(3{x), (3.69) 

where x{x) is an mi x mi matrix function. According to (3.1) and (3.67), 

respectively, f3 and 6 are boundedly differentiable on [0, /], and so x is also 
boundedly differentiable. Now, equalities (3.68), (3.69) and the first relations 
in (3.4), (3.58) yield that x' = (i.e., x is a constant). It follows from (3.1), 
(3.66) and (3.69) that x(0) = Imi, and therefore x = /^^, that is, 9 = (3. 
Thus, (3.59) follows directly from (3.66). ■ 

In a way which is quite similar to the proof of (3.57) and (3.58) we get 

y = -v*^, 7V = 0. (3.70) 

In order to recover 7 from /3, we recall that (5 is boundedly differentiable and 
/3/3* = /mj. Therefore, an 1712 x m matrix function 7, which is continuous, 
peacewise differentiable (with the differentiability dissappearing, possibly, on 
a finite set only), having bounded left and right derivatives everywhere on 
[0, /] and satisfying relations 

/37* = 0, 77* >0, 7(0) = 7(0), (3.71) 
can be easily constructed. Since 

/37* = /37* = 0, 77* = /m„ 77* >0, 7(0) = 7(0), (3.72) 
the matrix function 7 admits representation 

7 = X7, XX* >0, x(0) = /^„ (3.73) 

where x is continuous and peacewise differentiable. In view of (3.73), the 
second equality in (3.70) can be rewritten in the form 

= -x7'7*(77*)-\ x(0) = /^„ (3.74) 

which uniquely defines x. Hence, taking into account Corollary 3.12, Propo- 
sition 3.13 and formula (3.58), we see that we have obtained a procedure to 
solve the inverse problem. 
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Theorem 3.14 Let the potential v of the Dirac system (1.1) on [0, I] be 

bounded {i.e., let (2.25) hold). Then v can be uniquely recovered from the 
Weyl function in the following way. First, we recover $i using formula 
(3.53). Next, we use formula (3.59), where the operator Sx = PxSP* and 
S is given by (3.37), to recover /3. We recover {from /3) matrix functions 7 
and X via formulas (3.71) and (3.74), respectively. The matrix function 7 
satisfies the equality 7 = X7. Finally, we get v from the equality 

v{x) = ^'{x)-i{xy. (3.75) 

The last statement in this section is a Borg-Marchenko-type uniqueness the- 
orem, which follows from Theorems 3.11 and 3.14. More precisely, we will 
discuss the Weyl functions (/? and (p of two Dirac systems. Matrix functions 
associated with Ip are written with a "hat" (e.g., v, $1). 

Theorem 3.15 Letip and(p be Weyl functions of two Dirac systems on [0, /] 
with bounded potentials, denoted, respectively, by v and v. In other words, 
we assume that 

max(||v(x)||, \\v{x)\\) < M, 0<x <l. (3.76) 
Suppose that on some ray — cQz (c e M, 3^2; > 0) the equality 

\\ip{z) - ^{z)\\ = Oie^'"-^) {Qz 00) for all < r < Tq < Z (3.77) 
holds. Then we have 

v{x) = v{x), 0<x <ro. (3.78) 

Proof. Since, according to (2.15), Weyl functions arc non-cxpansive, we see 
that the matrix function e'"^^^^ {^Lp{z) — (p{z)) is bounded on the line '^z = 
M + e {e > 0) and the inequality 

lle-^''-^ {ip{z) - (p{z)) II < 2e2'"l^l {Qz > M > 0) (3.79) 

holds. Furthermore, formula (3.77) implies that e~^^^^ {^ip{z)—ip{z)) is bounded 
on the ray ^z — cQz. Therefore, applying the Phragmen-Lindelof Theorem 
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with the angles generated by the hne — M + e and the ray — cQz 
(^z > M + s), we get 

||e-2'"" {ip{z) - ip{z)) II < Ml for "iz > M + e. (3.80) 

Because of formula (3.38), of its analog for (p, $i and of the inequality (3.80), 
we have 

r 

e2'(^-'')"($i(x) - ^i{x))dx = Oilf^iz)) C^iz) ^ +oo) (3.81) 

uniformly with respect to ^{z). Hence, the left-hand side of (3.81) is bounded 
in C and tends to zero on some rays. Thus, we obtain 

f e2>(^-'^>($i(x) -$i(x))dx = 0, i.e., $i(x) = $i(x) (0<x<r). 
Jo 

(3.82) 

Since (3.82) holds for all r < ro, we get $i(x) = $i(x) for < x < tq. In 
view of Theorem 3.14, the last identity implies (3.78). ■ 

4 Inverse problem on the semi-axis and evo- 
lution of the Weyl function for NLS 

1. In this section, we again consider Dirac system (1.1) on the semi- axis 
[0, oo) with a bounded potential v: 

\\v{x)\\<M fora;eM+. (4.1) 

According to Definition 2.5, the Weyl function of system (1.1) such that (4.1) 
holds is given by the right-hand side of (2.5). Prom Proposition 2.2 we see 
that the Weyl function (p is unique. Moreover, it follows from the proof of 
Corollary 3.12 that 

/•oo 

ip{z) = 2iz e2-^$(a:)(ia:, e-^(^+^)-$(a:) G L^,,^,, (4.2) 

^0 



23 



and for any < Z < oo we have 

^{x) = {Q<x< I). (4.3) 

Since $i = we get that $i(x) does not depend on I ior I > x. Compare 
this with the proof of Proposition 4.1 in [15], where the fact that E{x, t) (and 
so $i) does not depend on / follows from the uniqueness of the factorizations 
of operators Sf'^. See also Section 3 in [4] on the uniqueness of the accelerant. 
Using (4.2), (4.3) and Theorem 3.14 we get our next result. 

Theorem 4.1 The unique Weyl function of Dirac system (1.1) satisfying 
(4.1) admits representation 

POD 

ip{z)^2iz e''^^^^{x)dx, e-'^^^,{x) e L^^,^^ {rj>M). (4.4) 
Jo 

The procedure to recover $i and v from tp is given in Theorem 3.14- 
2. The well-known "focusing" nonlinear Schrodinger (fNLS) equation 

'ivt + i{vx:c + '2vv*v) = (vt := ^u) , (4.5) 

where v{x,t) is an mi x m2 matrix function, is interesting from a mathe- 
matical point of view and plays an essential role (in its scalar, matrix and 
multicomponent forms) in numerous applications (see, e.g., [2,44]). The fNLS 
equation (4.5) is equivalent (see [11,48] and references therein) to the zero 
curvature equation 

Gt-F, + [G,F]^0, [G,F]:^GF-FG, (4.6) 

where the m x m {m = mi + 777.2) matrix functions G{x,t,z) and F{x,t,z) 
are given by the formulas 

G = izj + jV, F = i{z^j - izjV - (K + jV')/2) , (4.7) 

and j and V are defined in (1.2). The zero curvature representation (4.6) of 
the integrable nonlinear equations appeared soon after the seminal Lax pairs 
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(see a historical discussion in [11] and original papers [1,25,47]). Furthermore, 
condition (4.6) is the compatibility condition for the auxiliary linear systems 

y. = Gy, yt = Fy, (4.8) 

and in the case of the fNLS equation the first of those systems is (in view 
of (4.7)) the skew-self- adjoint Dirac system (1.1). The compatibility condi- 
tion was studied in a more rigorous way and the corresponding important 
factorization formula for fundamental solutions was introduced in [38,40]. 
It was proved in greater detail and under weaker conditions in [34]. More 
specifically, we have the following proposition. 

Proposition 4.2 [34]- Let m x m matrix functions G and F and their 
derivatives Gt and exist on the semi-strip 

V = {(x, i) : < X < oo, < i < a}, (4.9) 

let G, Gt and F be continuous with respect to x and t on T>, and let (4.6) 
hold. Then we have the equality 

u{x, t, z)R{t, z) = R{x, t, z)u{x, 0, z), R{t, z) := R{0, t, z), (4.10) 

where u{x, t, z) and R{x, t, z) are normalized fundamental solutions given, 
respectively, by: 

u^ = Gu, u{0,t,z) = lm; Rt = FR, R{x,0,z) = lm. (4.11) 
The equality (4.10) means that the matrix function 

y — u{x, t, z)R(t, z) = R{x, t, z)u{x, 0, z) 

satisfies both systems (4.8) in V. Moreover, the fundamental solution u ad- 
mits the factorization 

u{x, t, z) = R{x, t, z)u{x, 0, z)R{t, z)'^. (4.12) 

To derive the evolution Lp{t, z) of the Weyl functions of the skew- self-adjoint 
Dirac systems y-xix, t, z) = G{x, t, z)y{x, t, z), we rewrite (4.12) in the form 

u{x, t, z)-^ = R{t, z)u{x, 0, z)-^R{x, t, z)-^. (4.13) 
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Note that an additional parameter t appears now in the functions discussed 

in Sections 2 and 3 (see, e.g, u{x,t,z) and (fi{t,z)). Wc partition R{t,z) into 
blocks R — {Rij}^j^i (and Ru is an mi x mi block here). 

Theorem 4.3 Let an mi x m2 matrix function v{x,t) in T> he continuously 
differentiable together with its first derivatives and let v^x exist. Assume that 
V satisfies the fNLS equation (4.5) as well as the following inequalities: 

sup 11^(2;, t) II < M, sup ||t'a;(a;, t)|| < 00 for each < t < a. (4.14) 
{x,t)ev xeT&+ 

Then the evolution (fi{t, z) of the Weyl functions of Dirac systems yx — Gy, 
where G has the form (4.7), is given {for '^{z) > M) by the equality 

cpit,z) = {R2iit,z) + R22it,z)if{0,z)){Rnit,z) + Ri2it,z)if{0,z)y\ (4.15) 

Proof. Since V = V* and jV = —Vj, it follows from (4.7) that F{x, t,z)* — 
—F{x,t,z). Hence, using (4.7) and (4.11) we get 

^(i?(x, t, z)*jRix, t, z)) = Rix, t, zTijFix, t, z) - F{x, t, z)j)R{x, t, z) 
= R{x,t, z)* {i{z^ - z^)Im + {z + z)V{x,t) + iVx{x,t)j)R{x,t, z). (4.16) 

In view of (4.14) and (4.16), for each < t < a there are numbers Mi^2{t) 
(Ml > M > 0, M2 > 0) and the corresponding quarterplane 

Vi^{z: Qz> Ml, ?ftz > M2} C C, (4.17) 

such that the inequality 

^{R(x,t,z)*jR{x,t,z)) <0 (4.18) 

holds in Vi. Inequality (4.18) and the initial condition R{x, 0, z) = 1^ imply 
R{x, t, z)*jR{x, t, z) < j, or, equivalently, 

{R{x, t, zfY^jRix, t, z)-^ > J, ze V,. (4.19) 

Let V{z) satisfy (2.3) and thus satisfy (2.1). Let V{x,t,z) (sometimes we 
will also write V{z), omitting x and t) be determined by the equality 

V{x, t, z) := R{x, t, z)-^V{z). (4.20) 



26 



In view of (4.19) the matrix function V satisfies (2.1) in Vi. 

Because of the smoothness conditions on v, we see that G and F given by 
(4.7) satisfy the requirements of Proposition 4.2, that is, (4.13) holds. Using 
(4.13) and (4.20), we see that 

u{x,t,z)-^V{z) = R{t,z) 

where 

(f>{x,t,z)^[0 Im,]u{x,0,z)-'V{z)(^[lm, O] u{x,Q, z)-'V{z)'^ \ (4.22) 

Here, (2.9) yields inequalities 

det ( [I^, 0] u{x, 0, z)-'Viz)^ ^ 0, det ( [l^, O] u{x, t, z)-'V{z)) ^ 0. 

(4.23) 

According to (4.21) and (4.23) we have (in the quarterplane Vi) the equality 
[0 Im,]uix,t,z)-'Viz)(^[lm, 0]uix,t,z)-'V{z)') 

= (i?2i(t, z) + i?22(t, z)<P{x, t, z)) (i?ii(t, z) + i?i2(t, z)(t){x, t, z)) ~\ (4.24) 

Taking into account Definition 2.1, we see that the left-hand side of (4.24) 
belongs J^{x,t,z) and (f){x,t,z) G A/'(a;,0,z), where is defined in (4.22). 
Therefore, (4.15) follows from (2.5) and (4.24) (see also (2.24)), when x 
tends to infinity. Although, we first derived (4.15) only for T>i, we see that 
it holds everywhere in Cm via analyticity. ■ 

Remark 4.4 Theorems 4-i o-nd 4-3 can be applied to recover solutions of 
the fNLS. Theorems on the evolution of the Weyl functions also make up the 
first step in proofs of uniqueness and existence of the solutions of nonlinear 
equations via the ISpT method (see, for instance, [35]). 
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6{x,t,z)_ 



[V 0]u{x,0,z)-'V{z), (4.21) 
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